Isomorphs in the phase diagram of a model liquid without inverse power law repulsion 
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It is demonstrated by molecular dynamics simulations that liquids interacting via the Buckingham 
potential are strongly correlating, i.e., have regions of their phase diagram where constant-volume 
equilibrium fluctuations in the virial and potential energy are strongly correlated. A binary Buck- 
ingham liquid is cooled to a viscous phase and shown to have isomorphs, which are curves in the 
phase diagram along which structure and dynamics in appropriate units are invariant to a good 
approximation. To test this, the radial distribution function, and both the incoherent and coherent 
intermediate scattering function are calculated. The results are shown to reflect a hidden scale 
invariance; despite its exponential repulsion the Buckingham potential is well approximated by an 
inverse power-law plus a linear term in the region of the first peak of the radial distribution function. 
As a consequence the dynamics of the viscous Buckingham liquid is mimicked by a corresponding 
model with purely repulsive inverse-power-law interactions. The results presented here closely re- 
semble earlier results for Lennard- Jones type liquids, demonstrating that the existence of strong 
correlations and isomorphs does not depend critically on the mathematical form of the repulsion 
being an inverse power law. 



INTRODUCTION 

Recently a series of papers has been published concern- 
ing so-called strongly correlating liquids and their physi- 
cal properties [lHfl . Liquids that exhibit these strong cor- 
relations have simpler thermodynamic, structural, and 
dynamical properties than liquids in general. A strongly 
correlating liquid is identified by looking at the corre- 
lation coefficient of the equilibrium fluctuations of the 
potential energy U(rx, ...,r^-) and virial W{r\, ...,rjy) = 
— 1/3 Yli r i ' V r; [/(ri, ...,rjv) [8] at constant volume: 
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Here brackets denote averages in the NVT ensemble 
(fixed particle number, volume, and temperature), A de- 
notes the difference from the average. The virial W gives 
the configurational part of the pressure Q , 

P V = Nk B T( Pl , . . . , pat) + W(n, . . . , r N ). (2) 

Strongly correlating liquids are defined [l[ as liquids that 
have R > 0.9. 

The origin of strong WU correlations was investigated 
in detail in Refs. @,13] for systems interacting via the 
Lennard- Jones (LJ) potential: 

v{r) = 4e 



(3) 



The fluctuations of W and U are dominated by fluctu- 
ations of pair distances within the first neighbor shell, 
where the LJ potential is well approximated by an ex- 
tended power Law (eIPL), defined as an inverse power 
law (IPL) plus a linear term Q: 



The IPL term gives perfect UW correlations, whereas 
the linear term contributes little to the fluctuations at 
constant volume: when one pair distance increases, oth- 
ers decrease, keeping the contributions from the linear 
term almost constant (this cancellation is exact in one 
dimension). The consequence is that LJ systems inherit 
some of the scaling properties of the IPL potential - they 
have a "hidden scale invariance" 4,9]. Prominent among 
the properties of strongly correlating liquids is that they 
have "isomorphs" , i.e., curves in the phase diagram along 
which structure, dynamics, and some thermodynami- 
cal properties are invariant in appropriate units 0|. 
The physics of strongly correlating liquids was briefly re- 
viewed recently in Ref. 0- 

Since the LJ systems consists of two IPL terms, it is 
perhaps tempting to assume that a repulsive (inverse) 
power law is necessary for the hidden scale invariance de- 
scribed above. In the present paper we use the modified 
Buckingham (exp-six) pair potential to show that this is 
not the case. The Buckingham potential was first derived 
by Slater from first-principle calculations of the force be- 
tween helium atoms [l(|. Buckingham later used this 
form of the potential to calculate the equation of state 
for different noble gases [ll| . The Buckingham potential 
has an exponential repulsive term, while the attractive 
part is given by a power law (l2l . ITii ] : 



f(0 = e( s ^«p[a(l-^ r )]- 5 § s 



(5) 



Veip L {r) = Ar~ n +B + Cr. 



(4) 



Here e is the depth of the potential well and r m specifies 
the position of the potential minimum. The parameter a 
determines the shape of the potential well. The Bucking- 
ham potential is better able to reproduce experimental 
data of inert gasses than the LJ potential [14H16T ]. but 
is also computationally more expensive (unless look-up 
tables are utilized ffl). 
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FIG. 1. (a) The correlation coefficient, R, as a function of 
temperature on isocfiores for single-component Buckingham 
(SCB), single-component Lennard- Jones (SCLJ), and argon. 
For SCB and SCLJ argon values were used for all potential 
parameters QJ], and R was calculated directly from Eq. JTJ. 
For argon R was calculated from experimental data [7J|] as de- 
scribed in Refs. [J, 0]. The correlations are strongest for state 
points with both high density and high temperature, and the 
difference between the Buckingham and the LJ potential is 
small. The correlation coefficient R > 1 for low-temperature 
20.0 mol/L argon is of course unphysical and either caused 
by an uncertainty in the experimental data or the approxi- 
mations applied in the calculation of R (see Refs. [3. l20l| for 
details), (b) The value of 7 (Eq. ©) plotted versus tem- 
perature for the same systems as in (a). For argon 7 was 
calculated from experimental data [3 using Eq. (J7J). 7 de- 
creases slowly for increasing temperatures, except when the 
correlation coefficient is low (R < 0.9). 



All simulation data in this paper were obtained from 
molecular dynamics in the NVT ensemble. The samples 
contained 1000 particles. The simulations were set up 
by instant cooling from a high temperature state point 
followed by an equilibration period, to ensure the simu- 
lations were independent from each other. The simula- 
tions were performed with the RUMD molecular dynam- 
ics package [l7j |. which is optimized for doing computa- 
tions on state-of-the-art GPU hardware. 



ingham systems are very similar to those of argon and the 
LJ system (dotted line in Fig. [T]). This is a first indica- 
tion that the actual functional form of the repulsive part 
of the potential does not have to be an inverse power law 
in order for a system to exhibit strong WU correlations. 

Another interesting property of the fluctuations is the 
quantity 7 defined [4|, [5J as 



(AWAU) 
((AC/) 2 ) ' 



(6) 



When a system is strongly correlating (R is close to one), 
AW ss 7AC/. For IPL potentials 7 is constant and equal 
to n/3 and R = 1. For non-IPL potentials, however, 
R < 1 and 7 may chan ge w ith temperature and density 
as seen in Fig. [lib) flL 1 1 Ol) . Especially for R < 0.9, we 
find 7 changing rapidly. The curves are similar for the 
20.0 mol/L SCB and SCLJ systems, except for a vertical 
offset. 

Fluctuations in U and W are of course only directly 
accessible in simulations. For experimental systems one 
must revert to the use of thermodynamic quantities that 
reflect the fluctuations in U and W. For instance, the 
configurational part of the pressure coefficient (3y — 
(d(W/V) I dT)v and the configurational part of the iso- 
choric specific heat per unit volume dy — (d(U /V)/dT)y 
can be used to to calculate 7 for argon as follows [H, [H, Hj] : 



(7) 



The values of 7 for argon obtained in this way are plot- 
ted in Fig. [ljb) , and the agreement with the Buckingham 
systems is good. This confirms that the Buckingham po- 
tential produces more accurate predictions of experimen- 
tal argon data than the LJ potential. 

Interestingly, low density argon has a higher correla- 
tion coefficient than high density argon. This is the op- 
posite of what is found for the Buckingham and the LJ 
potentials. Furthermore, the buckingham data are in bet- 
ter agreement with the argon data at low density than at 
high density. At the present we do not have any expla- 
nation for these observations. 



CORRELATIONS IN SINGLE-COMPONENT 
BUCKINGHAM LIQUIDS 

To compare the simulations with experiments [l8j , ar- 
gon parameters from Ref. [14| were used; a — 14.0, 
r m = 0.3866 nm, e/fcs = 123.2 K. As can be seen in 
Fig-UJa), the single-component Buckingham (SCB) liq- 
uid is strongly correlating ( R > 0.9) in parts of the phase 
diagram, particularly at high densities and/or tempera- 
tures. The correlation coefficients (Eq. (JTJ) of the Buck- 



ISOMORPHS IN BINARY BUCKINGHAM 
MIXTURES 

Strongly correlating liquids are predicted to have iso- 
morphs, which are curves in the phase diagram along 
which structure, dynamics, and some thermodynamical 
properties are invariant in appropriate reduced units 0, 
6]. Introducing reduced coordinates as fj = p 1 / 3 ^, two 
state points (1) and (2) are defined to be isomorphic if 
pairs of microscopic configurations with same reduced co- 
ordinates (f^ = f( 2 ^) have proportional configurational 



Boltzmann weights: 

e -y(r< 1 \...,r«)/kBr 1 =Ci2e -U(r^,...,rP)/k B T 2 _ ^ 

Here the constant C\i depends only on the two state 
points and Eq. (jSJ) is required to hold to a good ap- 
proximation for all physically relevant configurations |6j . 
An isomorph is a curve in the state diagram for which 
all points are isomorphic (an isomorph is a mathemat- 
ical equivalence class of isomorphic state points). The 
isomorphic invariance of structure, dynamics, and some 
thermodynamical properties - all in reduced units - can 
be derived directly from Eq. (JSj) [5]. Only IPL liquids 
have exact isomorphs, but it has been shown that all 
strongly correlating liquids have isomorphs to a good ap- 
proximation (Appendix A of Ref. 0]). 

Among the thermodynamical properties that are iso- 
morphic invariant is the excess entropy, S ex = S — Sideai , 
where Sideai is the entropy of an ideal gas at the same 
temperature and density. In the following, isomorphic 
state points were generated by utilizing that the quantity 
7 in Eq. ^ can be used to change density and temper- 
ature while keeping the excess entropy constant 0, Q : 

By choosing the density of a new isomorphic state point 
close to the density of the previous isomorphic state 
point, the temperature of the new state point can be cal- 
culated from the fluctuations by combining Eq. ([5]) and 
Eq. (O 5] . In this way a set of isomorphic points can be 
obtained from one initial state point. 

The predicted isomorphic invariance of the dynamics 
is most striking in viscous liquids, where the dynam- 
ics in general depend strongly on temperature and den- 
sity. To demonstrate that a systems interacting via the 
Buckingham potential have isomorphs, we study what we 
term a Kob- Andersen binary Buckingham (KABB) mix- 
ture with potential parameters being the same as for the 
original Kob- Andersen binary L J (KABL J) mixture [2l[ : 
£aa = 1-0, r mi AA = tAB = 1-5, r. m%A B = 0.8^2, 
e B B = 0.5, r„^BB = 0.88-^2. A 4:1 mixture (A:B) was 
used with a = 14.5. The potentials were truncated and 
shifted at rff = 2.5r TO)ii /^2. 

One of the predicted invariants on an isomorph is 
the structure of the system. To test this prediction, 
the radial distribution function in reduced coordinates 
g(r) = g{p 1 ^r) was plotted for isomorphic state points 
(Fig. [H[a)). The structure is invariant for the large (A) 
particle pair correlation function to a very good approx- 
imation. For the AB and BB pairs the structure is less 
invariant. However, when a comparison is made with 
Fig. [2th), it is clear that g(r) for the AB and BB pairs is 
still more invariant on an isomorph than on an isotherm 
(note that the density variation on the isomorph is larger 
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FIG. 2. The radial distribution functions g(r) for simulations 
of the KABB mixture. Both graphs are in reduced units where 
f = p 1 ^ 3 r. (a) g(f) for isomorphic state points and the three 
different particle combinations. The structure is invariant on 
the isomorph for the AA particle pairs, but for the AB and BB 
pairs the structure is less invariant, (b) g(f) for isothermal 
state points of smaller density variation. The structure is not 
invariant on the isotherm for any of the particle pairs. 



than on the isotherm). This situation is similar to what 
is found for the KABL J system 

To investigate the dynamics of the systems, the inco- 
herent intermediate scattering function F s (q, t) is plotted 
in reduced units in Fig.[3ta) and Fig.^b). The presence 
of a plateau in F s shows that the system is in a viscous 
state, where the dynamics are highly state point depen- 
dent. The large difference in F s for the two isothermal 
state points confirms this (dashed lines). For the iso- 
morph all F s data collapse more or less onto the same 
curve, showing that the dynamics are indeed invariant 
to a good approximation on an isomorph. In contrast 
to the radial distribution functions, the invariance holds 
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FIG. 3. Incoherent intermediate scattering function for the 
A (a) and B (b) particles of the KABB system. The time 
is given in reduced units (i = p 1//3 T 1//2 t) and the q-vector 
is kept constant in reduced units: q A = 7.25(p/1.2) 1//3 and 
qB = 5.5(p/1.2) 1 ' 3 . The solid lines represent isomorphic state 
points, while dashed lines show isothermal density changes 
for comparison. The dynamics are to a good approximation 
invariant on an isomorph when expressed in reduced units, 
especially when compared to the isotherm. In contrast to 
g(f), this holds for both the A and the B particles. 



well for both types of particles. 

To investigate the invariance in dynamics further, the 
coherent intermediate scattering function was calculated 
(Fig. The coherent intermediate scattering function 
was calculated from the spatial transform of the number 
density p(q) 0. In order to obtain good results, it is 
necessary to average over time scales that are 10-15 times 
longer than what is usual for the intermediate scattering 
function. This is the reason that there are less state 
points shown for the coherent-, than for the incoherent 
intermediate scattering function. The data confirm that 
the dynamics are invariant on the isomorph, especially 
when compared to the isothermal density change (dashed 
lines). However, the invariance seems to hold slightly 
better for the AB and BB parts, which is the opposite of 
what is seen for the structural invariance. 

For systems described by a generalized LJ potential 
consisting of two IPL terms, the invariance of the struc- 
ture leads to a prediction for the shape of an isomorph 
when plotted in the U-W plane @ (generalized LJ po- 
tentials are a sum of inverse power laws). Since the re- 
pulsive term in the Buckingham potential is described 
by an exponential function, it is not possible to derive 
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FIG. 4. Coherent intermediate scattering function for dif- 
ferent particle pairs on an isomorph in reduced units (t = 
p!/3yi/2^ The solid lines represent isomorphic state points, 
while dotted lines show isothermal density changes for com- 
parison. Again, the q-vector is kept constant in reduced 
units: q AA = 7.34(p/1.2) 1/3 , q AB = 8.01(p/1.2) 1/3 and 
q BB — 6.01(p/1.2) 1//3 . Also F(q,t) is invariant on the iso- 
morph. Contrary to what is seen for g(r), the invariance 
holds better for the AB and BB parts. 



an exact equation that describes the isomorph in terms 
of U and W. Figure [SJa) shows that isomorphs for the 
KABB system agree well with the prediction for the 12- 
6 LJ system if a = 14.5 (this value of a was chosen to 
demonstrate this feature). For a — 13.0, there is a signif- 
icant difference with the predicted shape at higher den- 
sity and temperature. Figure shows the isomorphs 
for both values of a after scaling U and W by the same 
isomorph-dependent factor, demonstrating the existence 
of a master isomorph Q. This shows that master iso- 
morphs exist not only in generalized LJ systems where 
they can be justified from analytical arguments [6| . 
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FIG. 5. (a) Plot of the potential energy per particle versus 
virial per particle for the KABB system. The solid lines are 
the predictions of the isomorph shape for the 12-6 LJ poten- 
tial 0. For these predictions the initial state points with 
p — 1.2 were used as reference point. Since a new value 
of 7 was calculated for every state point, 7 is not constant 
on the isomorphs, but changes approximately 10% along the 
isomorphs. The shape of the KABB isomorph agrees very 
well with the predicted shape for the 12-6 LJ potential for 
a — 14.5. For a = 13.0, the shape is different, (b) The same 
data now scaled with Wq defined as the virial at U — [6J. 
The isomorphs scale onto each other, forming a so called mas- 
ter isomorph for each value of a. 



THE INVERSE-POWER-LAW (IPL) 
APPROXIMATION 



As mentioned in the introduction, a generic explana- 
tion for the existence of strong correlations and iso- 
morphs in non-IPL systems, is the fact that some pair 
potentials can be well approximated by an eIPL (Eq. [4]) 
as shown in Fig. [5] Putting this explanation to a test, it 
was recently demonstrated that structure and dynamics 
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FIG. 6. The figure shows how the Buckingham potential 
(a = 14.5) can be approximated by an extended IPL po- 
tential (eIPL). The red dotted line is the IPL approximation 
obtained using the parameters obtained below in Fig. [7] The 
difference of the IPL approximation and the Buckingham po- 
tential (dashed green line) is more or less linear in the first 
peak of g(r). By subtracting this linear term from the IPL 
term the eIPL approximation is found (dashed blue line). 



of the KABLJ system can be reproduced by a purely re- 
pulsive IPL system even in the viscous phase [22j ■ In the 
following we demonstrate that this procedure works also 
for the KABB system, despite its non-IPL repulsion. 

Following Pedersen et al. [22[ , we assume that the Kob- 
Andersen IPL (KABIPL) system used to approximate 
the KABB system has the form 



TPL 



(r) = Ae l3 ( f- 

1 i i 



(10) 



where the parameters eij and (Jij are the Kob- Andersen 
parameters for the different types of particles and the 
constants A and n are independent of particle type. 

For IPL liquids it is known that W = (n/3)U, so in 
principle the value of n could be calculated from 7 deter- 
mined from the WU fluctuations (Eq. ©). For non-IPL 
liquids however, there is a slight state point dependence 
of 7, so instead the slope of an isochore was used to de- 
termine n (Fig. E(a)) making use of the identity Q 



7 



tdW\ 
\~dUJv 



(11) 



For the Buckingham potential with a — 14.5 we obtained 
7 = 4.904 and n = 14.71. This is lower than the 7 = 5.16 
which was found for the 12-6 LJ potential [22j. This 
is also consistent with the data in Fig. [ljb) where the 
SCB system has a lower value of gamma than the SCLJ 
system. 
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FIG. 7. Three isochoric state points were used to obtain 
the parameters for the IPL potential, (a) The value of n 
was determined by linear regression to the mean values of 
the virial and the potential energy(marked by yellow crosses, 
(b) The method used to find the value of A of Eq. (TT2J. 
U IPL = X^>j e ij ( a i3/ r ij) n was calculated from configura- 
tions drawn from the KABB simulations and plotted against 
the energy obtained during the simulations; the value of A 
was then obtained from the slope of the mean energies (again 
marked by yellow crosses). 



From Eq. (fTU)) it follows that the total internal energy 
of the IPL system can be written as 



U 



IPL 



i>j 



(12) 



The scaling factor A was determined from the slope of the 
mean values of the energies in a U,U IPL plot (Fig. EZKb)), 
where U IPL is given by Eq. (fl"2"|) evaluated on configura- 
tions from simulations of the KABB mixture [22| ■ Using 
these parameters, simulations of the KABIPL systems 
were performed and the results were compared with the 
results of the KABB system. In Fig. [3] the incoherent 
intermediate scattering function of the two systems is 
plotted for comparison. The KABIPL reproduces the 
dynamics of the KABB system very well. It should how- 
ever be noted that in spite of the good reproduction of 
the dynamics, the KABIPL had a stronger tendency to 
crystallize than the KABB system at the two lowest tem- 
peratures due to the absence of attractive forces. The 
good agreement shown in Fig. [5] only holds if both sys- 
tems are in the same (supercooled) state. 

From the fluctuations in the potential energy one can 
calculate the excess isochoric specific heat using [8] : 



Gy =C V - ^Nk B = 



am 2 ) 

k B T 2 



(13) 



In Fig. HJa) Cy is plotted for different isochores cal- 
culated from KABB and KABIPL simulations. The 



FIG. 8. The incoherent intermediate scattering function 
F s {q,t), (q = 7.25) of the KABB and KABIPL simu- 
lations for isochoric state points with p = 1.2, T = 
0.42,0.44,0.46,0.50,0.6,1.0. The IPL potentials reproduce 
the dynamics of the Buckingham potential over a significant 
temperature range. 



heat capacities for the two systems follow each other 
closely, although there is a small and systematic differ- 
ence increasing with density. This is similar to what was 
found for the KABLJ system 22j, but the deviations 
are slightly larger for the KABB mixture. Figure Mjo) 
shows that the excess heat capacity to a good approx- 
imation obeys density scaling, Cy* /N = ,f(p 7 /T) and 
Rosenfeld-Tarazona scaling, C^/N = 5 (p)T" 2 / 5 @ - 
again in good agreement with results for the KABLJ sys- 
tem n . 



CONCLUSION 

The Buckingham potential has been shown to be 
strongly correlating like the Lennard- Jones potential. In 
spite of its exponential repulsion, the Buckingham poten- 
tial's dynamics and heat capacity can be closely approxi- 
mated by a purely repulsive IPL system. In particular the 
system has good isomorphs in the phase diagram. These 
findings are very similar to those found for Lennard- Jones 
systems. We conclude that the existence of strong corre- 
lations and isomorphs is not dependent on the repulsion 
being an inverse power-law. 



ACKNOWLEDGMENTS 

The centre for viscous liquid dynamics "Glass and 
Time" is sponsored by the Danish National Research 
Foundation (DNRF). 



7 



2.2 



u 



1.6 



1.0 




"B A 



-§- 



\ 



o KABB p = 1.128 
□ KABB p= 1.200 
a KABB p- 1.272 



KABIPLp = 1.12 8 



KABIPL p = 1 
KABIPLp = 1.2' 



2C0 



72 



12 3 4 
T 



0.5 



T/p' 



FIGi 9. (a) The ccnflgnrattonai part of tare intensive isoch|oric 

spec ific heat Cf? /N as a function of temperature. Three 



chores were simulated of the KABB and the KABIPL syste 
At low density the agreement between the two systems is fi 
good, but for higher densities the differences become 
(b) The same data plotted versus T/p 1 where 7 = 4.904 



collapyu on a single curve, wttcli hIiowh ttat ctensity t 



dat; 

ing worka. The function (lAZp 1 /T) 2/b waa fitted to the Jata 

(dashed line) , showing that Rosenfeld-Tarazona scaling is |also 
obeyed. 



[1] U. R. Pedersen, N. 
Pyre, Phys. Rev. 



P. Bailey, T. B. Schroder, J 
Lett. 100, 015701 (2008), 



10. 1103/PhysRevLett. 100.015701 



[2] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schre^ 
J. C. Pyre, J. Ch em. Phys. 129, 184507 (2008) 



10.1063/1.2982247 



[3] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schroder 



J. C. Pvre. J. Ch em. Phys. 129. 184508 (2008). 



10.1003/1.2982249 | 



doi 



[4] T. B. Schroder, N. P. Bailey, U. R. Pedersen, N. Gian 
J. C. Pyre, J. Ch em. Phys. 131, 234504 (2009) 
10.1063/1.3265955| 
[5] N. Gnan, T. B. Schroder, U. R. Pedersen, N. P 



ley, J. C. Pyre, J. Chem. Phys. 131, 234504 (2009), 



[6 



[7; 



[11 

[12; 

[13 



ms. 
rly 
lai ger. 
The 



m 

[15; 
fie 



c. 



doi 



doi 



TT7 

[18 

[19 

[20 

[21 

[22 



3oT 



[23 



10.1063/1.3265957 

T. B. Schroder, If. Gnan, U. R. Pedersen, N. P. Bai- 
lcv. J. C. Pvre. J.l chem. Phys. 134, 164505 (2011), |doi 



10.1063/1.3582900 



U. R. Pedersen, r| Gnan, N. P. Bailey, T. B. Schroder 
J. C Pyre, J. Nc n-Cryst. Solids 357, 320 (2011), |doi 



10.1063/1.3582900 



M. P. Allen, P 
liquids (Oxford University 
T. B. Schroder, V 
vaerd, J. C. Pyre. 



Phys. Rev E 80, 041502 (2009), |doi 
10.1103/PhysRev$.80.041502j 



J. C. Slater, Phys. Rev. 32, 349 (1928), doi 



10.1103/PhysRev. 



52.34£ 



R. A. BuckineiiL 
(1938), |doi:10 109i i 



i. McMahan. M. Ross. Phys Rev. B 



P. A. Young, A 
24, 5119 (1981) 
L. Koci, R. Ah 
J. Phys.: Condehs. Matter 19 , 016206 (2007), |doi 



coi:10.1103/PhysRevB. 24.5119 



lja, A. Belonoshko, B. Johans son. 



10.1088/0953-8984 



/19/1/016206 



E. A. Mason, W. I ). Rice, J. Chem. Phys. 22, 843 (1954), 



doi:10. 1063/1. 174C 



10.1103/PhysRev. 



M 



Roskilde TJniver, 
http : // rumd . org 
E. W. Lemmon 
■mophysical proper 
WebBoo k, NIST 
69, http : //webbo<|>k 
P. Cosl ovich. C. IV 



(2009), ioi:10.106. 



EPAPS document 
tailed information 
data., http: //www 
W. Ko b. H. C. _ 
(1994), |doi:10.110; 
U. R. Pedersen 
Phys. Rev. Lc 



10. 1103/PhysRevI ,ett . 105. 157801 



doi:10.1080/00268 179809483145 



3ai- 



doi 



Tildesley, Computer simulations of 

Press, Oxford, 1987) 
R. Pedersen, N. P. Bailey, S. Tox- 



Proc. R. Soc London A 168, 264 
/rspa. 1938.0171 



200 



J. E. Kilpatri ck, Vf. E. Keller, E. F. Ham mel, Phys. Rev 
97, 9 (1955), _ 
A. A. Abrahamscln 



doi: 1 0.1103/PhysRev.97.9 



Phyi 



s. Rev. 130, 693 (1963), doi 



L30.693 

~rty Molecular Dynamics package, 



O. McLinden, P. G. Friend, Ther- 
ies of fluid systems, NIST Chemistry 
Standard Reference Patabase Number 
.nist .gov 

Roland. J Chem. Phys. 130, 014508 



/1.3054635 



E-PRLTAO- 100-033802 contains de- 
on the analysis of experimental argon 
. aip . org/pubservs/epaps .html 
4-ndersen. Phys. Rev. L ett. 73, 1376 
/PhysRevLett. 73.1376| 

T. B. Schroder, J. C. P vre 
tt. 105, 157801 (2010), |doi 



Y. Rosenfeld, P. r ?arazona, Mo l. Phys. 95, 141 (1998) 



